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ABSTRACT 

A catalog of ionized gas velocity fields for a sample of 30 spiral and irregular galaxies of the 
Virgo cluster has been obtained by using three-dimensional optical data. The aim of this sur- 
vey is to study the influence of high density environments on the gaseous kinematics of local 
cluster galaxies. Observations of the Ha line by means of Fabry-Perot interferometry have 
been performed at the Canada-France-Hawaii, ESO 3.6m, Observatoire de Haute-Provence 
1.93m and Observatoire du mont Megantic telescopes at angular and spectral samplings from 
0"4 to 1"6 and 7 to 16 km s^^. A recently developed, automatic and adaptive spatial bin- 
ning technique is used to reach a nearly constant signal-to-noise ratio (S/N) over the whole 
field-of-view, allowing to keep a high spatial resolution in high S/N regions and extend the 
detection of signal in low S/N regions. This paper is part of a series and presents the inte- 
grated emission-line and velocity maps of the galaxies. Both Ha morphologies and kinematics 
exhibit signs of perturbations in the form of e.g. external filaments, inner and nuclear spiral- 
and ring-like structures, inner kinematical twists, kinematical decoupling of a nuclear spiral, 
streaming motions along spiral arms and misalignment between kinematical and photometric 
orientation axes. 

Key words: Galaxies: clusters: individual: Virgo - Galaxies: kinematics - Galaxies: veloc- 
ity field - Galaxies: dynamics - Instrumentation: Interferometry - Techniques: high angular 
resolution 



1 INTRODUCTION 

Dense environments such as the cores of galaxy clusters are known 
to give to disc galaxies photometric and spectral properties that may 
considerably differ from those of galaxies found in lower density 
environments, like the clusters periphery or the field. 

Manifestations of environmental effects in local galaxy clus- 
ters are illustrated by the morphology-density relation (Dressier 
1980; Binggeli, Tammann & Sandage 1987), the neutral hydro- 
gen deficiency of spiral galaxies in the cluster core (Chamaraux, 
Balkowski & Gerard 1980; Giovanelli & Haynes 1985; Gavazzi 
1987; Cayatte et al. 1990; Bravo-Alfaro et al. 2000) but also at 
large clustercentric distances (Solanes et al. 2001; but see San- 
chis et al. 2004) or the occurence of peculiar Ha morphologies 
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for spirals within the highest galaxy surface density regions (Moss 
& Whittle 2000; Koopmann & Kenney 2004; Vogt et al. 2004). 
Centres of nearby clusters are dominated by spheroidals, gas-poor, 
non star-forming discs (lenticulars) and ellipticals while lower den- 
sity regions are gas-rich, spiral discs dominated. Many works have 
furthermore brought evidence that the properties of cluster spirals 
strongly evolve with the redshift (see Poggianti 2004 and refer- 
ences therein). This can be seen in the evolution of the morphology- 
density relation (Dressier et al. 1997; Fasano et al. 2000; Smith et 
al. 2005) and of galaxy colours (Butcher & Oemler 1978, 1984). 
Distant, intermediate redshifts (0.2 ^ z < 0.8) clusters are pop- 
ulated by a larger number of star-forming spirals rather than by 
passive lenticulars. These results rise the question of a possible 
morphological transformation of spirals into SOs since a redshift 
of 2 ~ 0.5 and of the evolution of the star formation activity of 
galaxies (e.g. Couch & Sharpies 1987; Balogh et al. 1999; Dressier 
et al. 1999; Adami et al. 2000). An active star formation period has 
been rapidly quenched (e.g. Poggianti et al. 1999), leaving the star- 
forming galaxies into a poststarburst phase, a probably important 
step during their possible evolution towards early-type objects (e.g. 
Tran et al. 2003). The population of massive elliptical galaxies does 
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not vary as a function of redshift (Dressier et al. 1997; Postman et 
al. 2005; Smith et al. 2005) and their formation appears to go back 
to the earliest epochs > 2 — 3, e.g. Ellis et al. 1997; van Dokkum 
et al. 1998; van Dokkum & Franx 2001). 

The recent theoretical works armed with N-body or hydro- 
dynamical numerical simulations have furthermore considerably 
improved our understanding of the physical processes that could 
drive the evolutions of the stellar populations and the morphologi- 
cal properties of cluster galaxies. Among several theories, two pre- 
ferred invoked processes are tidal effects and ram pressure strip- 
ping (RPS). On one hand, simulations of tidal effects can account 
for a slow morphological change within a time-scale of few Gyrs 
when a galaxy enters a cluster and undergoes several gravitational 
interactions with other galaxies and with the cluster tidal field. In- 
stances of such models are the "galaxy harassment" (Moore at al. 
1996, 1998, 1999) or the "tidal heating" (Gnedin 2003a,b). In these 
models, low density galaxies such as low surface brightness discs 
lose almost all their stars; their remnants appear to become dwarf 
spheroidal systems and their stripped stars the diffuse intracluster 
light. Galaxy mergers could participate to the formation of central 
giant ellipticals and cD galaxies during the earliest phase of the 
cluster formation (Merritt 1985; Dubinski 1998). On another hand, 
ram pressure stripping (Gunn & Gott 1972) could explain the Hi 
deficiency of today spirals, and perhaps why the gas reservoir has 
been emptied and the star formation stopped in distant spirals. Here 
the hot gas of the intra-cluster medium (ICM) exerts a pressure on 
the gaseous disc of a galaxy, acting like a wind that can efficiently 
sweep it in several tens of Myrs (Abadi, Moore & Bower 1999; 
Quilis, Moore & Bower 2000; Schulz & Struck 2001; Vollmer et 
al. 2001 ; Roediger & Hensler 2005). Eccentric orbits allow galaxies 
to go deeper through the cluster core and thus to be more efficently 
stripped (Vollmer et al. 2001) and could explain the presence of Hi 
deficient spirals observed at large cluster-centric radius (Solanes 
et al. 2001). Because this mechanism does not involve gravitational 
forces, at a first order it cannot influence the (old) stellar population 
of a disc. At a second order, some simulations claim that the rise of 
gas surface density when a disc passes close to the cluster core and 
the gas shocks induced by the return of expelled gas into a parent 
galaxy can enhance the formation of stars (e.g. Vollmer et al. 2001). 
This is nonetheless probably not sufficient to explain the observed 
morphological segregation and the two above-mentionned mecha- 
nisms may be both at work to explain the observations (Poggianti 
et al. 1999). 

The kinematical and dynamical properties of cluster galaxies 
have also been the subject of several works although they have been 
less studied than their spectro-photometric ones. In local clusters, 
optical and Hi rotation curves are used to study the TuUy-Fisher 
relation of spirals (e.g. Schommer et al. 1993; Dale et al. 1997, 
1999; Verheijen 2001) or the shape of the curves and its connec- 
tion with the environment (e.g. Amram et al. 1993, 1996; Speran- 
dio et al. 1995; Adami et al. 1999; Rubin et al. 1999; Dale et al. 
2001). Rubin et al. (1999) used a large sample of optical rotation 
curves of Virgo cluster galaxies and found that the curves harbour- 
ing the most disturbed shapes have a distribution of systemic veloc- 
ities very well correlated with that of the Virgo ellipticals and are 
more and more frequent (relative to the regular rotation curves) as 
the cluster-centric radius from M87 (at the cluster core) increases. 
Both these results suggest highly eccentric orbits for spirals (see 
also Adami, Biviano & Mazure 1998). Moreover, the perturbed 
curves are not found preferentially in Hi deficient galaxies, sug- 
gesting that the two mechanisms responsible for the deficiency and 
the observed kinematical disturbances differ. In distant clusters, ob- 



servations of galaxies currently concentrate on the study of their 
Tully-Fisher relation (e.g. Ziegler et al. 2003; Bamford et al. 2005) 
and it is still too early to learn about any evolution as a function of 
redshift of their kinematical disturbances. 

In this work, we are interested in studying the kinematical and 
dynamical properties of Virgo spiral galaxies. For that purpose a 
sample of 30 galaxies obtained using optical data at high angular 
and spectral resolutions has been observed. Typical questions that 
are still to be answered and that motivate our study are : how do the 
properties of gaseous velocity fields, their asymmetries (induced 
by e.g. a warp, a lopsidedness, a bar or spiral pattern) and of dark 
matter haloes depend on such a dense environment ? Do morpho- 
logically disturbed spirals exhibit kinematical perturbations ? How 
does the kinematics correlate with the star formation efficiency ? 
Long-slit spectroscopy data can hardly answer these questions be- 
cause they only inform us about the kinematics along a few direc- 
tions (mainly the photometric major axis) and are not appropriate 
to investigate the full kinematical structure of spiral galaxies. The 
recent advent of a new generation of 3D optical spectroscopic in- 
struments offers very promising perpectives for the stellar and/or 
gaseous kinematics of large samples of nearby galaxies spanning 
large ranges of luminosity and morphology (de Zeeuw et al. 2002; 
Emsellem et al. 2004; Bershady et al. 2005; Daigle et al. 2005a; 
Hernandez et al. 2005a). These instruments are e.g. fiber-fed or 
lenslet array integral-field spectrographs (Bacon et al. 2001; Ber- 
shady et al. 2004) or Fabry-Perot (FP) spectro-imagers (Gach et 
al. 2002). We want to extend the available data with a sample of 
gaseous velocity fields of Virgo cluster spirals using the Fabry- 
Perot interferometry technique. 3D optical data are a good compro- 
mise for the study of environmental effects in clusters due to a bet- 
ter distribution and resolution of the gas kinematics than for radio 
data. Virgo cluster spirals are indeed Hi deficient, often showing Hi 
discs smaller than the optical ones (Cayatte et al. 1990). Moreover 
the neutral hydrogen often lack their central parts, regions where 
the molecular gas is mainly present (Sofue et al. 2003a,b). As an 
additional motivation to our study, the 3D observations of gas will 
be compared to results of N-body and/or hydrodynamical simula- 
tions, as already done for a few Virgo galaxies (see e.g. Vollmer et 
al. 2004). It could help to constrain the evolutionary scenario of pe- 
culiar objects in the cluster (e.g. ram pressure stripping, encounters 
of galaxies). 

This article is the second from a series. The first paper was 
dedicated to NGC 4438, a prototype of a galaxy interacting with 
the ICM and with companions (Chemin et al. 2005). This article 
presents a catalog of 2D Ha emission-line and kinematical maps 
for 30 Virgo galaxies observed with the Cigale and FaNTOmM^ 
instruments. The article presents the sample of the galaxies, the 
principles of the data reduction and of the adaptive binning in sec- 
tions |2| and |3| Sections |4] and |5| briefly present the results and a 
discussion. Section |6|resumes the results presented in the article. 
Comments for the individual galaxies are given in Appendix lAll 
and their kinematics is presented in Appendix lA2l 

The data will be soon available from a Fabry-Perot data-base. 
They can be obtained from the authors before upon request. 
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Table 1. Virgo galaxy targets parameters 
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Notes on columns: (1-3): Galaxy VCC, NGC-IC and Messier names. (4-5): Right ascension and declination (J2000). (6): Morphological type, from 
RC3 (de Vaucouleurs et al. 1991). (7): Total B-band magnitude, from RC3. (8): Systemic velocity (in km s'^), from LEDA (Paturel et al. 1991). 
(9): Projected distance to M87 (in degree). (10): Revised region of membership, as defined in Gavazzi et al. (1999): A = Cluster A, B = Cluster B, E = 
East cloud, M = M cloud, N = North cloud, S = Southern extension. (1 1): Hi deficiency parameter from Gavazzi et al. (2005). "The value is uncertain 
for NGC 4438 due to uncertainties on its morphology and optical diameter. 



2 DEFINITION OF THE SAMPLE 

The Virgo cluster of galaxies is the nearest cluster from the Milky 
Way (D ~ 16 Mpc) and is therefore very well appropriate to ob- 
serve in details the galaxy properties. The Virgo Cluster Catalog 
(hereafter VCC, Binggeli, Sandage & Tammann 1985; Binggeli, 
Popescu & Tammann 1993) counts ~ 1400 members with a domi- 
nating population of dwarf ellipticals. The cluster is not yet relaxed, 
as shown by the presence of irregular structures in either kinemati- 
cal or X-ray data (Binggeli et al. 1987; Bbhringer et al. 1994; Shi- 
bata et al. 2001). The 3D structure of the cluster is very particular 
(see e.g. Gavazzi et al. 1999), articulated around a main substruc- 
ture centered on Messier 87 (cluster A) and with several extensions. 

The selection criterium was to observe targets among the 
brightest spiral galaxies of the cluster. This turns out to select 
objects having a total magnitude in the B— band brighter than 

= 12. Other galaxies were observed because of their pre- 
sumably perturbed major-axis kinematics, as derived from long-slit 
spectroscopy (Rubin et al. 1999), their peculiar morphology (e.g. 
IC 3225) or their membership within a galaxy pair (e.g. IC 3583). A 
concern was to specifically observe galaxies in the inclination range 
between 25° and 80° (with the exception of NGC 4438) in order 
to avoid problems of disc opacity at optical wavelength for high 



inclinations and of determination of the circular velocities for low 
inclinations. The Fabry-Perot catalog is composed of 30 galaxies 
which locations are shown on Fig.|2|and properties given in Tab.0 
The morphological types of the targets are distributed through 3 
so/a, 1 Sa, 3 Sab, 3 Sb, 7 Sbc, 4 Sc, 3 Scd, 3 Sd, 1 Sdm, 1 Sm and 
1 Im galaxies. The targets are mainly located in the cluster core 
where environmental effects are supposed to be more important, 
but also in the cluster extension towards M49. Deep Ha imagery is 
available for all galaxies from this sample (Koopmann, Kenney & 
Young 2001; Boselli & Gavazzi 2002; Boselli et al. 2002; Gavazzi 
et al. 2002, 2003). 

Optical kinematical data for many individual objects or larger 
samples of Virgo galaxies have been presented in the past using 
long-slit spectroscopy of the Ha and [Nil] emission lines (e.g. 
Chincarini & de Souza 1985; Distefano et al. 1990; Sperandio et 
al. 1995; Rubin et al. 1999; Yoshida et al. 2004). Except for few 
individual galaxies (e.g. NGC 4254, Phookun, Vogel & Mundy 
1993; NGC 4321, Arsenault, Roy & Boulesteix 1990, Knapen et 
al. 1995, Canzian & Allen 1997; NGC 4522, Vollmer et al. 2000), 
the mapping of the 2D optical kinematics for several Virgo galax- 
ies has never been done before this study. The kinematics of the 
molecular and neutral atomic gas have nevertheless already been 
mapped at millimeter and centimeter wavelengths. The FP catalog 
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Table 2. Telescope and Instrument configuration 
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has 10 galaxies in common with the Virgo CO Survey catalog of 
(NGC 4212, 4254, 4501, 4535, 4536, 4548, 4569, 4579, 4654 and 
4689; Sofue et al. 2003a,b), 6 with the BIMA Survey Of Nearby 
Galaxies CO catalog (NGC 4321, 4450, 4535, 4548, 4569, 4579; 
Heifer et al. 2003), 15 with the Hi catalog (NGC 4237, 4254, 4321, 
4438, 4450, 4501, 4535, 4548, 4567, 4568, 4569, 4579, 4639, 4654, 
4689; Guhathakurta et al. 1988; Cayatte et al. 1990, 1994) and 
other works like the one and two dimensional Hi data presented 
in Warmels (1988a,b,c). All of these will allow a comparison of the 
kinematics of the different gaseous components. 

One has to note that not all galaxies of the sample lie at the 
same distance from us. For instance, accurate estimates based on 
the observations of Cepheid variable stars put NGC 4535 at a dis- 
tance of 16.0 ± 1.9 Mpc (Macri et al. 1999) and NGC 4639 at 25.5 
± 2.5 Mpc (Saha et al. 1997). The three-dimensional structure of 
the Virgo cluster has often been studied (Yasuda, Fukugita & Oka- 
mura 1997; Gavazzi et al. 1999; Fouque et al. 2001; Solanes et al. 
2002) but is neglected here and for simplicity reasons it is assumed 
a distance of 16 Mpc for the cluster and all its galaxies. Therefore, 
1" on the sky corresponds to ~ 78 parsecs at this distance. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Data acquisition 

The observations of the 30 galaxies were performed from 1999 to 
2005 at the ESO 3.6m (Chile), the 3.6m Canada-France-Hawaii 
(CFHT, Hawaii), the 1.93m Observatoire de Haute-Provence (OHP, 
France) and the 1.6m Observatoire du mont Megantic (OmM, 
Quebec, Canada) telescopes. Each of our visitor instruments, 
FaNTOmM (at ESO, CFHT and OmM) and Cigale (at OHP) is 
a package including a scanning Fabry-Perot interferometer, a se- 
ries of interference filters and an Image Photon Counting system 
(IPCS) based on a technology of GaAs amplifier tube having a 
high sensitvity and a quantum efficiency of ~ 28% (Gach et al. 
2002; Hernandez et al. 2003). Our focal reducer (named Cigale at 
ESO and OHP, MOS-FP at CFHT, Panoramix) is attached at the 
Cassegrain focus of the telescopes and is equipped with our visitor 
intrument. Table|2|and Table|3|give the observational set-up at each 
telescope and parameters of each galaxy (respectively). Among the 
30 galaxies, NGC 4654 is the only one that has not been observed 
with an IPCS but with the CCD detector of the MOS/FP module at 
CFHT. 

To reconstruct a spectral profile in a spatial pixel during an 
exposure, the spectral range (hereafter referred as FSR for Free 
Spectral Range) of a FP interferometer is scanned through a num- 
ber of channels (nfech)- This number must satisfy the Nyquist cri- 
teria, nfoch ^ 2 X F, where F is the interferometer finesse. The 
finesse is related to the resolving power R of the interferometer by 



R — Ps X F, where ps is the interference order at the scanning 
wavelength (As). The spectral range at As is given by FSR — ^ 
(in A) or FSR — (in km s~^), where c is the speed of the light. 
Columns 3,7 and 8 of Table|3]give the interference order (at Ha), 
the free spectral range (at As) and the number of channels of the 
observations. At the Ha wavelength, the finesses of the four used 
interferometers with interference orders of p = 609, 765, 793, 899 
and 1 162 are between ~10 and 23. This leads to resolving powers 
and spectral resolutions ranging between ~ 7950 or 0.83 A (p = 
793 interferometer) and 21000 or 0.31 A (p = 899 interferometer). 
As a comparison, the resolution of long-slit spectroscopy observa- 
tions (Rubin et al. 1999) is of the order of 1 A. 

During an integration, each channel is scanned many times 
within a typical time of 15 seconds. When the sky transparency is 
not excellent prior to an exposure, the integration time is set to 10 
seconds per channel in order to average more efficiently the effects 
of transparency changes. The important time for a FP observation 
is the total exposure time of a spectral sampling element, i.e. the to- 
tal integration time of one channel (toxp, column 9 of Table|3j. The 
average overall integration time of the whole sample is 3.8 min- 
utes per channel. The longest times are spent for galaxies observed 
at the smaller telescopes and for an object like NGC 4438, which 
has a very low surface brightness, diffuse and filamentary Ha mor- 
phology. To ensure an optimal wavelength calibration, a reference 
datacube of a Neon emission- line (A = 6598.95 A) is obtained 
before and after the observation of a galaxy. 



3.2 Data reduction 

The Fabry-Perot data were all reduced following a homogeneous 
procedure which proceeds as follows: integration and wavelength 
calibration of the raw datacube, spectral smoothing (banning fil- 
tering) of the wavelength-calibrated datacube, night-sky emission- 
lines (hereafter NSEL, generally OH lines) subtraction, treatment 
of the pixels signal-to-noise (S/N) and computation of the Ha in- 
tegrated emission and velocity maps. The A-calibration is done rel- 
atively to the above-mentionned Neon emission-line. Notice that 
no flux calibration of the datacubes is done. 

We refer to other papers for a more detailed description of the 
wavelength calibration (e.g. Laval et al. 1987; Amram 1991). The 
reduction package consists of IDL^ routines that were adapted and 
improved from the adhocw software (Boulesteix 1993). The major 
improvements are for the sky subtraction and the signal-to-noise 
ratio (S/N) treatment procedures which are explained in detail and 
illustrated in Daigle et al. (2005b). A brief summary of these major 
implementations for our data reduction follows. 

- Night-sky emission-lines subtraction: The often-used sky 
subtraction method of the adhocw software consists in modelling a 
single NSEL average profile which is extracted from f.-o.-v. areas 
free from galaxy emission and removed from the observation. 
Instead of this "classical" NSEL subtraction method, the new 
procedure constructs a NSEL datacube which is fitted to the whole 
field-of-view (f.-o.-v.) and then subtracted from the observed 
datacube. Modelling the sky emission for each pixel instead of a 
unique pixel model reduces the NSEL residuals and thus allows 
a better detection of galaxy emission in low-surface brightness 
regions (e.g. interarm, external regions...). Notice that if the 
sky-model fit cannot converge to an adequate solution, which can 
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Table 3. Observational parameters 



Galaxy 


Set-up 


Date 


P 


As 


Af 


FWHMf 


FSR 


nbch 


^cxp 


t 


Seeing 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




B2 


m 0^ 0^ 


/ yj 


uouy .uo 


6612 




380 71 


24 




120 






B2 


27-05-01 


793 


UJU i .\J 1 


6561 


12 


377 95 


24 


3 7 


88 8 






B2 




793 


(SSSl 7(^ 


6582 




379 14 


24 


5 2 


124 8 




Nnr 4.9^4. 


^ 


14-02-05 


765 


6615 47 


6621 


18 


395 03 


48 


5 


240 




Nnr 49Q4. 


CI 


07-04-03 


899 


(SS70 fi4 


6575 


20 


333 87 


48 


2 5 


120 




Nnr 49QS 


Y) 


Of. 04 09 


/ yj 






12 




24 


Z.J 






IC 3225 




04-04-02 


609 


6614 51 


6617 


^ I 


496.15 


40 


2 3 


92 




Nnr 4.^9 1 






899 


6597 17 


6605 


15 


335 22 


52 


5 


260 




Nnr 4^ SI 


B2 


27-04-03 


793 


UU 1 -J. J>J 


6614 




380 96 


24 


4 5 


108 




Nnr 440^ 


B2 


Z.O U*! — \JJl 


793 




6602 


12 


380 25 


24 


4.5 


108 




Nnr 44^8 


J) 


UU Wt uz. 


793 


6564 33 


6568 


12 


378 14 


24 


6.5 


156 




Nnr 44^0 


D 


07-04-02 


793 


UU(.7J . JU 


6607 


12 


380 51 


24 


2 5 


60 




Nnr 44^7 


B2 


10-03-03 


793 


(SS89 OQ 


6582 




379.16 


24 


5 2 


124 8 




Nnr 44Q8 


Bl 


\JJl \JJl VJU 


793 


6595 77 


6593 


10 


379 95 


24 


3 2 


76.8 




Nnr 4*^0 1 


Bl 


1 0-0^-00 


793 


6612 71 


6614 


1 1 


380 93 


24 


4.8 


115.2 




Nnr 4"^ 1 Q 


CI 




899 


6589 49 


6595 


18 


334 83 


48 


2 


96.0 




Nnr 4's^9 


Y) 


0^-04-09 

U J U'T uz. 


793 


UUL/U.OZ. 


6607 


12 


380 59 


24 


2.0 


48 




NGC 4535 


A 


06-03-03 


899 


6605.71 


6617 


15 


335.66 


52 


3.0 


156.0 




NGC 4536 


A 


14-03-04 


765 


6602.36 


6595 


18 


394.25 


48 


3.4 


163.2 


* 


NGC 4548 


A 


07-03-03 


899 


6573.42 


6585 


16 


334.01 


52 


2.1 


109.2 


* 


NGC 4567 


Bl 


12-04-99 


793 


6612.56 


6614 


11 


380.92 


24 


5.0 


120.0 


* 


NGC 4568 


Bl 


12-04-99 


793 


6612.56 


6614 


11 


380.92 


24 


5.0 


120.0 


* 


IC 3583 


D 


05-04-02 


793 


6587.32 


6588 


11 


379.46 


24 


5.0 


120.0 




NGC 4569 


A 


11-03-02 


765 


6557.64 


6569 


15 


391.58 


40 


3.8 


152.0 




NGC 4571 


Bl 


08-03-00 


793 


6570.27 


6574 


12 


378.48 


24 


5.0 


120.0 


* 


NGC 4579 


D 


04-04-02 


609 


6596.03 


6598 


10 


494.76 


40 


2.3 


92.0 




NGC 4580 


B2 


08-03-03 


793 


6585.42 


6584 


12 


379.35 


24 


5.5 


132.0 




NGC 4639 


D 


05-04-02 


793 


6584.89 


6588 


11 


379.32 


24 


3.3 


79.2 




NGC 4654 


C2 


13-02-02 


1162 


6585.48 


6588 


11 


258.89 


24 


2.5 


60.0 




NGC 4689 


A 


03-03-03 


899 


6598.16 


6595 


18 


335.27 


52 


3.2 


166.4 


* 



Notes on columns: (1): Observing configuration : see Tablel2lfor details on the telescope and instrument used during the observa- 
tions. (2): Date of observation (dd-mm-yy). (3): Interference order (at the Ha wavelength). (4): Scanning wavelength (A). (5-6): 
Interference filter central wavelength and FWHM (in A, at the ambient temperature). (7): Free Spectral range (at the scanning 
wavelength, in km s^^). (8): Number of channels. (9): Total exposure time per channel (in min). (10): Total scanning exposure 
time (in min). (11): is for a seeing < I'.'S, * for a seeing > 2'.'0 and < 4'.'0. 



happen when there are two competing NSELs or when the galaxy 
emission covers a very large majority of the f.-o.-v., the "classical" 
method is used. This has been done e.g. for NGC 4438 and NGC 
4579. 

- Datacube signal-to-noise treatment: An adaptive spatial bin- 
ning is used in order to increase the S/N of the data, instead of 
the often-used gaussian smoothing. For that purpose IDL routines 
were adapted to our Fabry-Perot data from existing routines de- 
veloped by Cappellari & Copin (2003) for the SAURON data (de 
Zeeuw et al. 2002). The method is a Voronoi tessellation and ba- 
sically consists in spatially binning pixels with their neighbouring 
pixels until a bin with sufficient S/N is obtained. Consequently, 
in regions of initial high S/N (Hll regions in the spiral arms, in 
the disc, in the galactic centre, ...), no spatial binning (or a mini- 
mum binning) is done and the angular resolution is kept maximal 
in these regions, contrary to the gaussian smoothing method. On 
another hand, because the binning occurs for pixels having a. S/N 
lower than the target S/N value (low surface brightness and dif- 
fuse gas in interam or outermost regions, ...), it allows to recover 
the spectral information over a larger part of a f.-o.-v. than the gaus- 
sian smoothing method. The data presented here have a minimum 
S/N of 5 over the field-of-view. An example of a binned velocity 



field is shown for the galaxy NGC 4254 in Fig. Q with spectra of 
bins (having different sizes) selected at some locations inside and 
outside the disc. 

For 3 galaxies (NGC 4405, NGC 4438 and NGC 4569), a 
gaussian filtering of 3 x 3 pixels FWHM (3x3 and 6x6 pixels for 
NGC 4438) was used instead of the Voronoi tessellation because 
the solution of this new procedure was not optimal for them. It can 
happen when the binning is done on a NSEL residual (and so not 
on the Ha line of the galaxy) in extremely low surface brightness 
regions (e.g. filaments of NGC 4438, external arm of NGC 4569). 

The CCD data of NGC 4654 had two additional pre- 
processing procedures. Indeed, the CCD data require a cosmic ray 
hits subtraction and a photometric correction, procedures which are 
not yet implemented in our IDL routines. Because each channel is 
scanned only one time during a CCD exposure with an integration 
of several minutes per channel, significant sky transparency varia- 
tions can occur during an exposure and a photometric correction to 
the data is needed to take into account the sky flux changes. This 
problem does not exist with a photon counting system because each 
channel is quickly scanned several times during an exposure so that 
the transparency variations are averaged from the beginning to the 
end of the overall integration. For the CCD observation of NGC 
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Figure 1. Ho binned velocity field of NGC 4254 obtained by the Voronoi 
tessellation. The contours represent the Ha integrated emission of the 
galaxy. For each spectra, a dashed line corresponds to the systemic velocity 
of NGC 4254. 

4654, the cosmic ray hits were removed and a photometric correc- 
tion was applied to the raw datacube as explained in Boulesteix 
(1993) before the wavelength calibration. 

The final steps of the data reduction consists in integrating the 
flux under the Hq emission line, deriving its barycentre, dispersion 
and underlying continuum level following the method described 
in Daigle et al. (2005b). We are finally left with an integrated Ha 
emission-line map and a Ha velocity field. All radial velocities pre- 
sented in this article are given in the heliocentric rest frame and are 
not absolute values. Notice that no systemic velocities are given in 
this article; their values which do not present any particular inter- 
est for this study are generally in very good agreement with the 
referred ones (column 8 of Tab. within less than 15 km s~^. 

Regions of galaxies that exhibit several emission lines in their 
profiles have their velocity fixed to the one of the highest intensity 
component (e.g. NGC 4438, Chemin et al. 2005). Mutiple spec- 
tral components are often found in interacting systems (Amram et 
al. 2004; Rampazzo et al. 2005). The FP continuum map roughly 
traces the underlying stellar continuum of the Ha line. It is not as 
deep or accurate as broad i?-band images because it is integrated 
over ~ 5 — 10 angstroms. The FP continuum images have actu- 
ally very little interest. For this reason, no continuum images are 
proposed in the article. Their use only concerns the recovery of the 
data astrometry (see next section). 

3.3 Data astrometry 

No reference to the World Coordinate System (WCS) have been 
obtained during the observations. It is important to know how a 
datacube is oriented in order to derive the major-axis position an- 
gle (hereafter P. A.) of a disc and the astrometric solution has to 
be found for the FP data. The koords task in the KARMA package 
(Gooch 1996) is used to recover it, using a i?-band map and a FP 
underlying continuum map as reference and target images respec- 
tively. The task provides a reference pixel, its angular size on the 
sky plane and a possible rotation angle to align the f.-o.-v. with the 



WCS axes. All the FP images presented in ADpendix lA2l have their 
vertical axis aligned with the North. The original FP f.-o.-v. are very 
well oriented with the WCS. The orientation angle of the FP f.-o.-v. 
vertical axis to the North is lower than 0.5°. The worst misalign- 
ments reach ~ 2° and are for 9 galaxies observed at OmM and 
OHP. The values of position angles of major axis given in Tab. |4| 
have been corrected of such strong misalignments. 



4 RESULTS 

4.1 Ha emission line and kinematical maps 

Figure|2|gives a global representation of all 30 velocity fields from 
this sample positioned at the location of galaxies in the Virgo clus- 
ter. Appendix lA 1 [ presents a brief description of the observed kine- 
matics for the galaxies of the sample. Appendix IA2| presents their 
Ha emission line and kinematical maps. 

The catalog is arranged with increasing order of right as- 
cension. For each galaxy, four images are displayed : a B-band 
image (Top-Left), a H-hanA image (Top-Right), the FP Ha inte- 
grated map (Bottom-Left) and the corresponding Ha velocity field 
(Bottom-Right). The B- and i/-band images are taken from the 
GOLDMine archive (Gavazzi et al. 2003). For IC 3225, the T^-band 
image is obtained from the 2MASS archive (Jarrett et al. 2000). All 
four maps are WCS oriented and have the same field-of-view. The 
scales of the 4 images is chosen so that both the large scale stel- 
lar morphology and Ha distribution and kinematics are displayed 
with great detail, with the noticeable exception of NGC 4438. The 
position-velocity (PV) diagrams are also displayed, using the de- 
rived rotation curve projected onto the sky plane (solid line). The 
datacube slice has a width of 3-4 pixels. No PV diagrams are pre- 
sented for NGC 4438 and NGC 4548 for the reasons exposed in 
A constant value of P.A. is adopted for NGC 4457, NGC 
4501 and NGC 4579 despite the fact their kinematical P.A. ap- 
peared to vary with radius (see i|4.3> . 

Preliminary results for IC 3225, IC 3583, NGC 4298, NGC 
4450, NGC 4532, NGC 4569, NGC 4579 and NGC 4639 have been 
presented in Chemin (2003) and detailed results for NGC 4438 in 
Chemin et al. (2005). This catalog presents the first optical kine- 
matical data ever obtained for the galaxies IC 3225, IC 3583, NGC 
4405 and NGC 4457. Forthcoming papers from this series will de- 
tail more quantitatively results on the whole sample and on individ- 
ual galaxies. 

The FP emission line maps are in agreement with the deep 
Ha images presented in Koopmann et al. (2001) and in the GOLD- 
Mine database (Gavazzi et al. 2003). The discs have morphologi- 
cal perturbations such as nuclear spiral (NGC 4579), nuclear, in- 
ner and outer pseudo-ring (NGC 4321, NGC 4639), off-plane dif- 
fuse filaments, arc-like structure and giant Hll regions (NGC 4438, 
NGC 4532, NGC 4569), highly truncated Ha discs (NGC 4405, 
NGC 4569, NGC 4580), lopsidedness (e.g. IC 3225, NGC 4298, 
NGC 4654), m — 2 modes such as a bar or two spiral arms, more 
complex spiral structure (NGC 4254, NGC 4535, NGC 4654). The 
signatures of kinematical perturbations are observed in the veloc- 
ity fields in the form of streaming motions along spiral arms or 
pseudo-ring (e.g. NGC 4254, NGC 4294, NGC 4321, NGC 4639, 
NGC 4654, NGC 4689), steep velocity gradients along a bar or a 
nuclear structure (e.g. NGC 4519, NGC 4536, NGC 4579), Z-shape 
of the velocities due to the presence of a bar (e.g. NGC 45 19, NGC 
4532, NGC 4654), decoupled nuclear spiral (NGC 4579), twist of 
the kinematical major axis (e.g. NGC 4501) or misalignment be- 
tween photometric and kinematical major axes (NGC 4254, NGC 
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Figure 2. A global view of the Virgo cluster galaxy velocity fields map. For clarity reasons, the angular scale of each galaxy has been enlarged by a factor of 
~ 9.5 w.r.t. the real scale. The locations of the 3 ellipticals M49, M86 and M87 are displayed using broad-band optical images. The colour scale is blue-darker 
shades (red-lighter shades) for the approaching (receding respectively) side of galaxies. See Appendix IA2l for the velocity scale of each object. 
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4519), solid-body rise of the velocity in truncated discs. Another 
(more subtle) perturbation is a kinematical lopsidedness that is only 
revealed by an in-depth study of the velocity fields with titled-ring 
models and a harmonic decomposition (see Fig.|5|and section ls!2l . 
All these perturbations are probably responsible for the observed 
irregularities of long-slit rotation curves (Rubin et al. 1999). 

4.2 Kinematical parameters fitting 

An important step of the pre-analysis is to determine the location 
of the kinematical centre, the orientation parameters of a galactic 
disc (its inclination i and angle of the major axis P. A.) and the 
recession velocity Vsys of the disc from a kinematical map. Except 
for the inclination that is only a scale factor of the rotation curve and 
iigya, all the parameters are supposed known before an observation 
in long-slit spectroscopy studies. The kinematical and photometric 
major axes are indeed supposed to be exactly the same and have 
a constant orientation across the galaxy, allowing the alignment of 
the slit with the photometric major axis of a disc in order to derive 
its rotation curve. The dynamical centre is furthermore supposed 
to be at the location of the photometric one. It will be shown in 
the next section that these arguments sometimes do not hold and 
can lead to significant velocity (and thus mass) discrepancies with 
respect to the reality. 

There are several ways to derive a rotation curve Kot from a 
velocity field. Some 2D-kinematics optical studies fit by eye the pa- 
rameters to the velocity field with efforts to lower the visual differ- 
ence between the rotation curves of the receding and approaching 
halves of a disc (e.g. Garrido et al. 2002). They also do not permit 
any variation of the parameters as a function of radius. 

In this work, it is decided to use a more robust, rigourous 
and automatic fitting routine that fully uses the power of two di- 
mensional maps and with a minimum of human intervention. We 
use the rotcur task (Begeman 1987, 1989) of the gipsy package 
(van der Hulst et al. 1992). It considers purely axisymmetric mo- 
tions by fitting the formula Hobs = Wsys + Kot- sin(i). cos(S) + 
Kad- sin(j). sin(0) to the observed velocity fields. Kot and Vrad 
are the circular and radial velocities in the galaxy plane. The terms 
i, Wobs and v^y^ refer to the azimuthal angle in the galaxy plane, 
the disc inclination, the observed l-.o.-s. and systemic velocities 
(respectively). The radial velocity is very frequently not fitted in a 
velocity field (e.g. Rogstad, Lockart & Wright 1974; Bosma 1981; 
Begeman 1989; Carignan & Beaulieu 1989; van der Hulst et al. 
1993; de Blok, McGaugh & van der Hulst 1996; Verheijen 2001) 
because it is considered very negligible with respect to the rota- 
tional term. Using a sample of 8 spiral galaxies observed in CO and 
Hi, Wong, Blitz & Bosma (2004) estimate that radial motions of in- 
flow or outflow only reach an upper limit of 3-5% of the circular 
velocity in only 3 velocity fields from their sample. It is therefore 
chosen here to neglect the Kad term during the fittings. The pa- 
rameters to be adjusted as a function of radius are the coordinates 
of the centre of mass, v^ys, i and the position angle P. A. of the 
major axis. This latter is the counterclockwise angle that the semi- 
major axis of the receding half does with the North axis. A good 
sampling of a ring implies its width to be set to a value which is at 
least the average seeing of the observation and > 2 times the pixel 
size. It leads to ring width of > 2" (for data collected at ESO and 
CFHT), > 3" (at OHP) and ~ 5" (at Megantic). Previously, the use 
of rotcur on Hq velocity fields was not very successful because of 
the poor spatial coverage limited to the brightest Hll regions. But 
as can be seen in Appendix I A2I the spatial coverage of these new 
data allows rotcur to converge to accurate results, as with Hi data. 



It is mainly due to the new generation of IPCS detector used and to 
the new adaptive binning technique that allows to get information 
of the diffuse Ha component in interarm regions (see i|3.2t . 

The sequence of the computation is similar to that usually used 
for Hi velocity fields (e.g. Carignan & Puche 1990; Verheijen & 
Sancisi 2001) and proceeds as follows: 

• Average values of Vsys and of the dynamical center coordi- 
nates are found from their fitted radial profile. As an initial guess, 
the photometric parameters are used. The coordinates of the pho- 
tometric centre were fitted to near-infrared (NIR) images obtained 
from the 2MASS and GOLDMine databases. The average seeing 
of the NIR images is l'.'6 and the accuracy on the photometric co- 
ordinates is less than 1". The position of the photometric centre is 
displayed with a cross in the figures of Appendix IA2I The photo- 
metric position angles of the major axis and the inclination were 
taken in previous optical kinematical studies (Rubin et al. 1999), 
in the LEDA database and the Uppsala General Catalog of Galax- 
ies (Nilson 1973). Notice that in this first iteration the inclination 
is kept fixed at the photometric value. Keeping the value of P. A. 
fixed or free does not influence the determination of Vgys and of the 
coordinates of the rotation centre. 

• Then, another least-square fitting is done with these new val- 
ues of dynamical centre coordinates and Vsys as fixed parameters 
and with i, P. A. and the rotation velocity as free parameters. An 
average value of i and P. A. is computed from their fitted radial 
profile. 

• A last fitting is done by having i and P. A. fixed by their new 
constant value, with the rotation velocity as the only free parameter. 

A residual velocity field is finally obtained by subtracting a model 
to the observation. This map is used to determine when an optimal 
set of parameters is reached, i.e. when the global average and the 
dispersion (r.m.s.) of the residual velocities are minimum. 

All points inside a half-sector of 30° (sometimes up to 50° 
for large inclinations) around the minor axis were rejected and a 
I cos(^)| weight was applied to the l.-o.-s. velocities during the fit- 
tings (Begeman 1989). Both allow to minimize the contamination 
of the closest points to the minor axis where important projections 
effets occur and where the circular velocity term is fitted with dif- 
ficulty (cos(6) 0). The average parameters are computed with 
some restrictions. First, one must avoid the inner regions where the 
velocities are often affected by non-circular motions. This is notice- 
ably the case for barred galaxies (Hernandez et al. 2005a). Their 
presence is clearly visible in the PV diagrams of Appendix IA2I 
where the rotation curves rarely fit the emission-line maxima in the 
central parts. Then, one should not take into account rings which 
contain too few points. Results for such rings can give unrealistic or 
highly uncertain results when compared with other well-populated 
rings. This basically corresponds to rings having a number of points 
of less than ~ 15% that of the most populated ring, i.e. for the very 
inner and outer rings (see Fig. |3| top panel). Finally, one should 
avoid the regions where only one half of the disc is observed, or 
where the number of points in the fitting is highly dominated by 
that of only one side of a disc. 

The two "Siamese" galaxies NGC 4567 and NGC 4568 ob- 
viously superpose along the line-of-sight. The fittings should be 
partly affected by pixels of overlap in their receding halves. The 
method to decide which of the pixels belong to NGC 4567 or NGC 
4568 was to fit models until the residual velocities are minimum in 
the region of overlap. As a result, it appears that the region of over- 
lap has minimum residual velocities when they are attributed to the 
NGC 4568 velocity field. This is fully consistent with the photom- 
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Table 4. Virgo galaxy targets parameters 



Galaxy 


P.A.p 


ip 


P 




ik 


Offset 




(1) 


(2) 




(3) 


(4) 


(5) 


(6) 




9fi5 


4S 


251 


±6 




4.3 


335 




255 


J J 


259 


±2 


4Q -1- fi 


1.3 


101 




108 


51 


106 


±2 


53 di 3 


2.0 


156 


Nnr 4.754. 


45 


28 


69 


±3 


31 di 6 


5.2 


406 


Nnr 47 Q4. 


335 


71 


335 


±3 


68 di 6 


4.1 


320 




320 


58 


316 


±3 


58 di 3 


1.6 


125 


IC 3225 


213 


75 


215 


±4 


(70 di 12) 


2.7 


211 


Nnr 4.17 1 


153 


32 


151 


±3 


38 di 7 


6.6 


515 




80 


49 


73 


±5 




2.6 


203 




200 


54 


199 


±3 


57 it 6 


2.2 


172 




355 


43 


351 


±7 


(49 i 17) 


3.6 


281 




82 


34 


88 ± 22 


(27 i 19) 


6.1 


476 




133 


59 


136 


±4 


63 it 4 


7.7 


601 




140 


59 


141 


±2 


58 it 6 


4.3 


335 


Nnr 4.5 1 Q 


325 


40 


357 


±4 


40 it 5 


7.3 


569 


NGC 4532 


340 


64 


355 


±6 


(63 it 14) 


10.3 


803 


NGC 4535 


180 


46 


180 


±2 


38 it 10 


2.2 


172 


NGC 4536 


310 


67 


301 


±2 


68 it 3 


9.6 


749 


NGC 4567 


85 


49 


90 


± 1 


53 it 4 


2.7 


211 


NGC 4568 


23 


67 


27 


±3 


67 it 3 


2.2 


172 


IC 3583 


180 


72 


178 


±5 


(77 it 11) 


4.4 


343 


NGC 4569 


23.5 


65 


23 


±4 


67 it 5 


4.0 


312 


NGC 4571 


220 


28 


220 


±4 


(37 it 17) 


3.2 


252 


NGC 4579 


95 


38 


91 


±3 


44 it 8 


1.8 


140 


NGC 4580 


165 


40 


157 


±3 


46 it 10 


0.9 


70 


NGC 4639 


304 


49 


311 


±2 


48 it 6 


2.4 


187 


NGC 4654 


128 


57 


125 


±2 


62 it 7 


2.4 


187 


NGC 4689 


163 


36 


165 


±2 


(36 it 14) 


2.1 


164 



Notes on columns: (1-2): Photometric position angle P.A.p and inclina- 
tion ip, (in degree). Taken from Rubin et al. (1999), LEDA, UGC. (3-4): 
Kinematical position angle P.A.^ and inclination i]^, (in degree). Num- 
bers into parentheses indicate that the kinematical inclination has not been 
chosen to obtain the rotation curve but the photometric value. (5-6): Sky 
projected offset between the photometric and kinematical centres (in arc- 
sec and parsec respectively). 



etry where NGC 4568 seems to be projected at the foreground of 
the pair. The parameters of both galaxies are thus accurately deter- 
mined. 



4.3 Comparison with pliotometric parameters 

Figure|3|illustrates examples of resuUs from the tilted-ring model 
fittings to the velocity fields. Table |4| gives the kinematical orien- 
tation parameters. The error-bar given in Tab. |4|is the r.m.s. (la- 
dispersion) of the radial profiles of the inclination or the position 
angle. Figure |4| shows a comparison between the photometric and 
kinematical inclinations and position angles of the major axis. A 
similar comparison is made for the position of the photometric and 
kinematical centres in Fig. |5| No fitting has been done for the ex- 
tremely perturbed galaxy NGC 4438 (as explained in Chemin et al. 
2005) and for NGC 4548 (lack of independent points because of 
too large bins). No position-velocity diagrams are shown for these 
two galaxies. 

Position angle of the major axis.- The position angle of the 
major axis is generally the best determined parameter, with a typi- 
cal dispersion of 3°. The largest dispersions are for very perturbed 
galaxies (IC 3583, NGC 4532) or objects having a peculiar kine- 
matics (NGC 4450). A comparison of the kinematical P.A. with 
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Figure 3. (Top) Radial profiles of kinematical inclination for NGC 4535 
(Left) and NGC 4568 (Right). A dashed horizontal line represents the mean 
inclination and two dashed vertical lines delineate the radial range within 
which the mean value is computed. The top panels display the number of 
points of each annulus. (Bottom-Left) Radial profile of kinematical position 
angle for NGC 4501. Blue (red) triangles are for the approaching (respec- 
tively receding) half of the disc and filled diamonds for both sides fitted 
simultaneously. (Bottom-right) Radial profile of the X- and Y-coordinates 
(open and full circles) of the kinematical centre for NGC 4298. 



the photometric one is done in Figure |4] and shows that there is a 
general good agreement between the two angles. A few galaxies 
have a misalignment of more than 10° with the photometry (e.g. 
NGC 4189, NGC 45 19). An analysis of their residual velocity maps 
confirm that the P.A. derived from their velocity field is the angle 
that describes better the gas motions in the disc than when using 
the photometric P.A. (see 35. 1> . 

For differences lower than 10° (23 galaxies among 28 when 
removing NGC 4457 which has a strong kinematical twist), the 
mean absolute difference between the two angles is 3°, which is 
exactly the typical error-bar on the kinematical P.A.. For those 
galaxies, the velocity curves derived from long-slit data should give 
a good representation of the rotational pattern, providing that the 
position of the photometric centre is also in agreement with the dy- 
namical one. 

A few velocity fields exhibit significant kinematical twists 
(NGC 4457, NGC 4501 and NGC 4579). Kinematical twists are of- 
ten associated with a warp outside the optical disc and thus with a 
change of inclination of the Hi disc (e.g. Sicotte & Carignan 1997). 
For NGC 4501 and NGC 4579, the twists occur in the galaxy cen- 
tres. The NGC 4579 twist is clearly linked to its nuclear spiral (see 
Aopendix lAlt which orientation differs by ~ 90° from that of the 
large-scale spiral arms. An example is shown in Fig. |3| (bottom) 
with NGC 4501. Such a twist was also observed in CO data (Wong, 
Blitz & Bosma 2004) and was attributed to an unkown, hidden bar. 
Its presence remains however to be confirmed by photometry. For 
those galaxies, a mean value is derived within a radial range where 
the angle variations are the smallest. It is this mean value which is 
given in Tab.|4]and used in the making of their PV-diagram. Notice 
that for NGC 4457, the given P.A. dispersion remains very large 
because the P.A. continuoulsy decreases with radius. For the two 
remaining galaxies, the P.A. estimates outside that central region 
of variation are in good agreement with the photometric values. 

Inclination.- The kinematical inclination is generally higher 
than the photometric value but remains in agreement with the pho- 
tometry within the errors. A large dispersion for the inclination il- 
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Figure 4. Comparison between photometric and kinematical orientation pa- 
rameters: absolute value of relative position angles versus the kinematical 
position angle (Left). Inclination (Right). 

lustrates the difficulty to accurately determine it from high resolu- 
tion velocity fields with a tilted-ring model, mainly due to the cou- 
pling of the circular velocity and the inclination during the least- 
square fitting (Kot- sin(z) term of the line-of-sight velocity). The 
largest dispersions are mainly for inclinations lower than ~ 50°, 
as already shown in Begeman (1989) with Hi data. Moreover, as 
for the position angle, the inclination of mildly perturbed discs are 
less accurate than those having a more regular morphology. Fig- 
ure |3| shows examples of inclination fittings for the galaxies NGC 
4535 (cr = 10°) and NGC 4568 (a = 3°). The kinematical inch- 
nations for IC 3225, IC 3583, NGC 4189, NGC 4351, NGC 4450, 
NGC 4457, NGC 4532, NGC 4571 and NGC 4689 are very uncer- 
tain (cr > 10°) or very inconsistent with the photometric value (see 
Tab.|4}. This is not surprising since all these galaxies are either mor- 
phologically perturbed discs or have a presumably low photometric 
inclination. The photometric value was thus used for the extraction 
of their rotation curve and the modelling of their velocity field. 

It is worth noting that m = 2 perturbations of the potential 
of a disc have been shown to also affect the results of inclination 
and P.A. fittings in tilted-ring models (Schoenmakers, Franx & de 
Zeeuw 1997; Schoenmakers 1999; Fridman et al. 2005). The pa- 
rameter is observed to oscillate around its mean physical value and 
it results to an increase of the errors on it, as probably observed in 
many galaxies of this sample. Notice finally that Hi velocity maps 
with lower resolution than the Ha observations logically give lower 
dispersion of the inclination (e.g. Guhathakurta et al. 1988; Wong 
et al. 2004), which parameters compare well with our FP estimates. 

Position of the rotation centre.- An example of a radial pro- 
file of the kinematical coordinates is displayed in Figure|3|for the 
galaxy NGC 4298. The coordinates show a large variation as a 
function of radius. This implies a large dispersion around the mean 
coordinates. In view of this quite large variation of the coordinates, 
getting the most accurate rotation curve with an optical velocity 
field should perhaps require to leave the centre vary as a function 
of radius. This is indeed shown to give the most symmetric rotation 
curves (Beauvais & Bothun 2001). The impact of the oscillations 
of the coordinates of the dynamical centre on a rotation curve is 
nevertheless not the purpose of this article. It will be studied in a 
forthcoming paper. As usually done, it has been decided here to 
keep the centre fixed and at a constant position. Keeping it constant 
is precisely what allows to further detect a kinematical lopsided- 
ness in the harmonic analysis that could reflect an environmental 
effect on a galaxy (see ^5.2> . 

Figure |5| shows the difference between the coordinates of the pho- 
tometric and kinematical centres (see also in Tab. |4| and the im- 
ages of Appendix IA2t . The [0,0] location corresponds to that of 



the photometric centre and a circle to its uncertainty. Its radius is 
thus given by the seeing of the NIR plates (l'.'6). Though there are 
very few points inside this circle, the positions of the two centres 
seem consistent with each other due to the large scatter often ob- 
served around the kinematical values (as seen with NGC 4298). 
The largest departures are for perturbed galaxies (NGC 4457, NGC 
4498, NGC 4519, NGC 4532) but also for large, barred galaxies 
(NGC 4321, NGC 4536). For all these objects, it is probably more 
difficult to determine the position of the photometric centre than 
for other more regular galaxies. Also, the extinction by dust near a 
galaxy nucleus could play a role in the determination of the photo- 
metric centre, making it harder to find though NIR images are used 
here for that purpose. 

Consequences of centering effects on a rotation curve have been 
discussed in Beauvais & Bothun (2001). We qualitatively find the 
same result as them that when choosing the photometric centre in- 
stead of the kinematical centre, the velocities of the very innermost 
annuli of a rotation curve can change drastically. We illustrate this 
point in Figure |5| (right panel) with the galaxies NGC 4536, NGC 
4567 and NGC 4689. These objects were chosen because the de- 
parture between the centres exceeds the NIR plates seeing of l'.'6. 
The expression ^^'''"^ displayed, where Vphot is the rota- 

tion curve derived by fixing the centre at the photometric position 
(using the kinematical inclination and P. A.). A systemic velocity 
corresponding to this centre was fitted and a corresponding residual 
velocity field was built. The r.m.s. of this residual map is generally 
by very few km s~^ greater than that of the residual map built us- 
ing Kot and the kinematical centre. The optical centre is thus not 
adequately the centre of rotation, though not being necessarily so 
far from it within the estimated error-bars. 

The graph of Figure |5| shows that the difference between the two 
curves always decreases as a function of radius. The curves only 
differ by ~ 2% in the outer regions. The flat part of a rotation curve 
is therefore not affected by a minor change of centre. In the in- 
nermost regions of barred galaxies, differences of up to 70-80% 
are observed. This difference is not significant because the rota- 
tion velocity is highly uncertain, whatever the choice of a centre. 
This is very likely due to important non-circular motions of the gas 
kinematics in the central parts of barred galaxies. For non-barred 
late-type spirals, the velocity of the innermost annuli differ by ~ 
5 to 25%. This value has to be compared with the formal error-bar 
measured by rotcur at these radii, which is less than 5% of Kot- 
Therefore, contrary to the barred objects, the difference is critical 
because the velocities at these radii are essential in the determina- 
tion of the disc mass-to-light ratio in mass models. The choice of 
the centre is important for this type of galaxies. 



5 A SHORT ANALYSIS 

This section briefly describes some goals that will be achieved with 
the high-resolution 3D optical data in future articles from this se- 
ries. 



5.1 Comparison with other Ha works 

As we wish to compare our derived rotation curves with results 
from earlier optical studies (Chincarini & de Souza 1985; Speran- 
dio et al. 1995; Rubin et al. 1999), it is first needed to ensure 
whether the FP measurements agree with the long-slit data, at 
least along the photometric major axis. An example is illustrated 
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Figure 5. (Left) Comparison between the kinematical and photometric co- 
ordinates. The horizontal and vertical bars are the rm.s. of the kinematical 
coordinates (resp.) and the circle centered on [0,0] is the mean seeing of the 
NIR plates, i.e. the error on the position of the photometric centre. (Right) 
Comparison between the rotation curves Vrot and Vpi^^t for NGC 4567, 
NGC 4536 and NGC 4689. These curves are respectively obtained by fix- 
ing the rotation centre at the kinematical and photometric centres. 

in Fig. |51 with the galaxy NGC 4189. Both long-sMt and FP rota- 
tion curves extracted along the major axis with a constant P. A. 
of 265° agree very well. However the average kinematical P. A. 
(251°) clearly differs from the photometric one. As a consequence, 
the rotation curve of the galaxy derived with the tilted-ring model 
differs from the long-slit one (bottom panel of Fig. |6}. To check 
that our result is a better representation of the kinematical pattern 
than the long-slit one, it is displayed on Fig.0the residual veloc- 
ities obtained by subtracting from the observation two models. A 
first velocity model is based on the tilted ring fitting to the FP ve- 
locity field (left panel) and another model is based on the smooth 
long-slit rotation curve with a P. A. value fixed to the photometric 
one (right panel). Notice that this smooth model does not take into 
account the last two points of the long-slit rotation curve and that 
it does not pass through the zero velocity. This does not change the 
following result. The figure shows residual velocities more homo- 
geneously distributed in the disc for our model, with a dispersion of 
16 km s~^ in the residuals, whereas there is an excess of residual 
velocities on both sides of the major axis for the long-slit model, 
with a worst dispersion of 28 km s~^ in the residuals. This signa- 
ture is the one expected in a residual map when the model uses an 
erroneous P. A. (see Warner, Wright & Baldwin 1973). 

As explained in i|4.3l the photometry and the kinematics are 
in a relative good agreement, implying that long-slit data should 
only occasionally give incorrect results, like with NGC 4189 when 
a strong misalignment is observed. This claim will be tested within 
a future work for all the Virgo sample. 

5.2 Non-axisymmetry of gaseous velocity fields 

Galaxies in clusters are sensitive to the environmental effects. Their 
structures (gas, stars, dark matter halo) can be perturbed by the 
cluster tidal field, gravitational encounters with other galaxies, 
galaxy mergers, accretion of gas, etc... Such external events can 
trigger internal perturbations in a stellar disc like bars or oval dis- 
torsions (e.g. Noguchi 1987; Boumaud & Combes 2002), spirals, 
warps (Huang & Carlberg 1997), polar rings (Bournaud & Combes 
2003) or lopsidedness (Boumaud et al. 2005). Ram pressure strip- 
ping is furthermore efficient in creating lopsided gaseous discs and 
off-plane annular-like gas structures (e.g. Schulz & Struck 2001; 
Vollmer et al. 2001). All these perturbations are likely responsi- 
ble for the formation of the kinematical asymmetries observed in 
our Ha velocity fields (Appendixes I A H and ! A2t and in other Hi or 
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Figure 6. Comparison between the rotation curves of NGC 4189. (Top 
panel) Long-slit rotation curve (Rubin et al. 1999, filled squares), smooth 
model of long-slit data (solid line) and rotation curve from the FP velocity 
field extracted along the photometric axis of 265° in a 2" width pseudo- 
slit (open circles). (Bottom panel) Rotation curve of NGC 4189 derived 
from the tilted-ring model fitting to the velocity field (filled squares and 
triangles) for a value position angle of the major axis of 25 1 ° . Blue (red) 
triangles are for the approaching (respectively receding) half of the disc. 
The solid line is the smooth model of long-slit data. 

Ha data (Swaters et al. 1999; Beauvais & Bothun 1999, 2001; Ver- 
heijen & Sancisi 2001), as well as in their optical rotation curves 
(Rubin et al. 1999; Dale et al. 2001). 

Franx, van Gorkom & de Zeeuw (1994), Schoenmakers et al. 
(1997) or Wong et al. (2004) made an analysis of Hi and CO ve- 
locity fields by means of Fourier coefficients in order to study the 
non-axisymmetry of the gas kinematics. A similar approach has 
been used with stellar kinematics of early-type galaxies by means 
of kinemetry (Copin et al. 2001; Krajnovic et al. 2005). The tech- 
nique decomposes the line-of-sight velocity along each ring of the 
tilted-ring model into harmonic terms Cm and s„i of order m fol- 
lowing 

Vobs = Co + ^ c„i cos(m\l/) + Sm sin(m>I>) 

where ^ corresponds to the angle in the plane of the ring, the coef- 
ficient Co to the systemic velocity of the ring, the first order term ci 
to more or less the rotation curve and all other terms to non-circular 
motions (axisymmetric and non-axisymmetric). Schoenmakers et 
al. (1997) showed that a galaxy perturbed by a distortion of the 
potential of harmonic order m exhibits terms of order m ~ 1 and 
771 + 1 in its velocity field. 

This technique has also been applied to ionized gas veloc- 
ity fields (Krajnovic et al. 2005; Simon et al. 2005) and reveals 
itself here to be successful in detecting the signatures of non- 
axisymmetry of the 2D Ha kinematics, as seen in Figure |8| This 
diagram displays the Fourier coefficients as a function of radius af- 
ter an harmonic expansion up to an order of 777 = 3 of the velocity 
field of NGC 4298, using fixed and constant orientation parameters 
(given in Tab.|4j, kinematical centre position and systemic velocity. 
The si term is not null and the S3 term appears to slightly oscillate 
around zero. There is no bar nor warp in NGC 4298 so that these 
variations could be due to the presence of a weak 777 = 2 spiral 
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Figure 7. Residual velocity maps of NGC 4189. (Left) Using the kinematical P. A. of 251° and the tilted-ring model rotation curve (Fig.|6]bottom). (Right) 
Using the photometric P. A. of 275° and the smooth model of the long-slit data (Fig.|6]top). The velocity range and the colour-scale are the same for the 2 
maps. 
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Figure 8. Harmonic decomposition of the velocity field of NGC 4298. 
Squares (triangles) correspond to the cosine (sine, resp.) terms. 



pattern. The most striking feature is the significant variations of the 
m = 2 kinematical terms. This can be caused by m = 1 or m = 3 
modes. This latter hypothesis can be rejected because the Ha mor- 
phology of NGC 4298 does not show signatures of m = 3 struc- 
tures like e.g. 3 spiral arms. There is no direct evidence of a "one 
spiral arm" structure in NGC 4298 so that the m = 1 mode could 
be caused by a kinematic lopsidedness. The presence of the galaxy 
NGC 4302, companion of NGC 4298 in the sky plane but also in 
radial velocity, could be at the origin of the observed asymmetries. 
Beauvais & Bothun (1999) also observed drifting centres in their 
FP observations. While a large number of galaxies appear to ex- 
hibit a morphological lopsidedness, as revelead by galaxy imagery 
(Boumaud et al. 2005) or Hi data (S waters et al. 1999; Noorder- 
meer et al. 2005), evidence of kinematical lopsidedness such as the 
one probably detected here in NGC 4298 seem to be rare in the 
WHISP sample of Noordermeer et al. (2005). 

A forthcoming work will apply the Fourier analysis to the 
Virgo cluster sample (as well as velocity maps of field galaxies) 
and will investigate how the kinematical perturbances can be re- 
lated to the cluster environment. 



6 CONCLUSION 

This article has presented Fabry-Perot interferometry observations 
of the ionized gas for a sample of 30 spiral and irregular galaxies 
of the Virgo cluster. This is the first time that a sample of Ha inte- 
grated emission-line maps and velocity fields are collected in order 
to study the environmental effects on the two-dimensional, high- 
resolution kinematics of cluster galaxies. Two new improvements 
in the reduction of the FP data consist in a better sky subtraction 
and an adaptive binning. Both these procedures allow to achieve an 



optimal signal-to-noise in the field-of-view and a better detection 
of the diffuse ionized gas in interarm or external regions. 

The galaxies exhibit morphological irregularities like nuclear 
spiral, nuclear and inner pseudo-ring, off-plane diffuse filaments 
and arc-like structure, truncated disc, lopsidedness, bar and spiral 
arms. The velocity fields exhibit typical kinematical perturbations 
like streaming motions along spiral arms, twist of the major axis, Z- 
shape of velocities due to the presence of a bar, decoupled nuclear 
spiral or misalignment between photmetric and kinematical major 
axes. Future papers from this series will present the rotation curves 
and mass models of the galaxies and investigate whether the asym- 
metries observed in the velocity fields have a relation with such a 
dense environment. 
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APPENDIX A: THE KINEMATIC AL Ha CATALOG OF 
VIRGO GALAXIES 

Al Comments for individual galaxies 

NGC 4189: The velocity field of this barred multi-armed galaxy 
is slightly perturbed. The position- velocity diagram reveals that the 
Ha peak in the nucleus has an almost constant velocity, which 
signature was already observed in the long-slit data (Rubin et al. 
1999). The kinematical and photometric major axes are not aligned. 

NGC 4212: The galaxy presents a mildly perturbed velocity field 
in the northern and western regions. The optical and near-infrared 
images clearly reveal an important dust obscuration effect in the 
galaxy. This may disturb the gaseous kinematics of NGC 4212, as 
particularly seen in the western region. The long-slit data showed 
velocities down to ~ —240 km s~^ (Rubin et al. 1999), which 
are not observed here. As in the CO data (Sofue et al. 2003a), the 
core exhibits a steeper velocity gradient than for the rest of the disc. 
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This could be due to the presence of a hidden bar. A careful analysis 
of the NIR image by an isophotal ellispses fitting indeed seems to 
reveal the presence of an elongated inner structure (see also Buta et 
al. 2005). 

NGC 4237: The velocity field shows no real signs of perturba- 
tions, except for a clump in the receding half to the South (R.A. 
~ 12''17'"12.53^ Dec. ~ +15°18'58"), which exhibits a lower 
velocity (~ 900 km s~^) than its surroimdings (~ 950 km s~^). 
Streamings of gas are also observed. 

NGC 4254: The prominent spiral structure is asymmetric in this 
well known-galaxy, as well as its Ha velocity field. Streaming mo- 
tions are observed along the spiral arms. Our observation of NGC 
4254 is in good agreement with the Ha and Hi velocity fields 
presented in Phookun et al. (1993). Notice a region at R.A. --^ 
12''18"'57.56", Dec. - +14°26'55", Uobs - 2535 km s"') hav- 
ing a higher velocity than its surroundings (~ 2500 km s~^). CO 
data show a steep velocity rise in the innermost arcseconds (80 km 
s~^). It could be due to a small bar detected by Kranz et al. (2001). 

NGC 4294: With NGC 4299, NGC 4294 belongs to an appar- 
ent pair. Streaming motions are seen along the spiral arms. The 
receding half appears slightly pertubed, with Hii regions having 
lower radial velocities than the local diffuse gas (e.g. at R.A. ~ 
12''21"16.92",Dec. ~+Il°3l'0l"). 

NGC 4298: This apparently flocculent galaxy (Elmegreen et al. 
1999) has a companion NGC 4302 only distant by 2.'3 (or 10.7 kpc 
projected) and with nearly the same systemic velocity (v^i'^^^ = 
1149 km s~^, Binggeli et al. 1985). According to Koopmann & 
Kenney (2004), the truncated gas disc is due to ram pressure strip- 
ping and the stellar asymmetry to a tidal interaction with a compan- 
ion. The velocity field appears however mildly perturbed. Notice 
that the southern arm presents velocities more redshifted than the 
nearby parts of the disc at the same azimuthal angles. This perhaps 
indicates streaming motions along it or, less probably, a locally 
warped arm. A steep velocity rise is detected in the PV-diagram 
in the core of the galaxy. 

IC 3225: This magellanic-type galaxy has an enhanced star for- 
mation on the NE, approaching half of its disc whereas very little 
ionized gas is detected on the SW half. The i3— band distribution 
shows an elongated stellar disc to the SW whereas the NIR im- 
age has a truncated disc in this direction. As seen in the position- 
velocity diagram, the velocity smoothly increases with radius. The 
velocity field thus mainly displays a solid body rotating pattern. 
This is typical of highly inclined optically thick discs, which pro- 
jection effects prevents from observing the true rotation (Goad & 
Roberts 1981; Giovanelli & Haynes 2002). 

NGC 4321: The well-known galaxy NGC 4321 (Ml 00) has of- 
ten been mapped with Fabry-Perot interferometry (Arsenault et al. 
1990; Knapen et al. 1995; Canzian & Allen 1997; Hernandez et al. 
2005b). Evidence for three different pattern speeds (nuclear struc- 
ture, large-scale bar, spiral) has been shown using the Tremaine- 
Weinberg method on this rather regular velocity field (Hernandez 
et al. 2005b), which also show signs of streaming motions (Knapen 
et al. 2000). The core of the galaxy has a nuclear ring-like structure. 

NGC 4351: The galaxy is slowly rotating. Its has a peculiar 
gaseous morphology, showing an arm-like structure to the North. 
The stellar distribution appears lopsided, also with an arm-like 
structure in the northern and western regions of the galaxy. The 
stellar asymmetry could result from a past tidal interaction, or even 
a merger. This galaxy is also a very good candidate for an ongoing 
RPS (Koopmann & Kenney 2004). 

NGC 4405: The Ha distribution is truncated and the velocity 
map displays a solid-body rotating pattern. This is tmderstandable 



in a framework where only the central regions lying in the rising 
part of a rotation curve are observed. 

NGC 4438: Notice the off-plane filamentary morphology of the 
ionized gas to the West and South of the galaxy disc, which was 
already seen in Kenney et al. (1995) and Kenney & Yale (2002). 
No Hi was detected in the disc (Cayatte et al. 1990). The ionized 
gas is the best tracer of the complex morphology and kinematics of 
the galaxy, which is presented in detail in Chemin et al. (2005). 

NGC 4450: Like in the Hi data of Cayatte et al. (1990), this ane- 
mic galaxy has a patchy Ha distribution. The Ha spiral appears 
highly wound around the centre. The velocity field appears signif- 
icantly dispersed in the galaxy inner parts as well as in the outer 
ones but to a lesser extent. The nucleus displays a steep gradient 
200 km s^'^), which result was already seen in long slit data 
(Rubin et al. 1999). Little CO gas is observed in NGC 4450 (Heifer 
et al. 2003). 

NGC 4457: The distribution and the kinematics of this galaxy are 

very peculiar. There is little ionized gas in the receding side, which 
is traced by only few star formation regions while the approaching 
side consists in a single Ha armlet. The core region surprisingly 
exhibits the most blue-shifted velocities observed in the galaxy. 

NGC 4498: The emission is mostly located in Hii regions along 
the stellar bar-like structure. The kinematics of NGC 4498 is asym- 
metric: the approaching half is typical of a constant velocity rota- 
tional pattern while the receding side has solid-body rotation. 

NGC 4501: The Ha velocity map of NGC 4501 (M88) appears 
regular with the exception of non-circular motions associated to 
the spiral structure, as particularly seen along the minor axis. NGC 
4501 exhibits one of the largest velocity range among all the galax- 
ies presented here, indicating a very massive object. Notice that 
part of the NW Ha spiral arm is missing in the approaching half 
of the disc, which is probably due to a dust obscuration effect and 
to a low transmission of the interference filter. The CO kinematics 
in the ±4" region around the photometric centre is complex, as re- 
vealed by double-peaks in the molecular gas emission (Onodera et 
al. 2004). The Ha PV-diagram also shows a complex central kine- 
matics. Cayatte et al. (1990) showed compressed contours of Hi 
distribution to the SE. 

NGC 4519: The galaxy morphology is perturbed, as well as the 
velocity field where streaming motions are revelead along the mul- 
tiple spiral arms. Notice the strong Z-shape of the velocities due to 
the presence of the bar. This galaxy has a misalignment between 
photometric and kinematical major axes. 

NGC 4532: The velocity map is very perturbed. It displays a 
strong Z-shape kinematical pattern, which has already been ob- 
served in the Hi velocity field (Hoffman et al. 1999). A systemic 
velocity component 2010 km s~^) is observed along the north- 
western edge of the Ha distribution. NGC 4532 has the highest Ha 
equivalent width in the sample of Koopmann & Kenney (2004). It 
is not surprising while regarding its perturbed stellar and gaseous 
morphologies, with prominent star forming regions (like e.g. the 
off-plane one at R.A. - 12''34'"21.02^ Dec. - -i-06°27'38") and 
stellar extensions. All the perturbations may result from a tidal in- 
teraction and/or gas accretion (Hoffman et al. 1999). Notice that an 
optical reflection of the galaxy was masked, giving artificial sharp 
edges to the shape of the galaxy in its eastern parts. 

NGC 4535: This barred galaxy has a low inclination and an ir- 
regular morphology, as mostly seen in the S-band and Ha images, 
with several spiral-like structures arising to the North and South in 
addition to the main spiral pattern. The velocity field is perturbed 
in the northern Ha loop and in the southern arm-like extension. 
Streaming motions are also observed along the arms. The PV dia- 
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gram shows a steep velocity rise associated with the bright Ha core. 
The central regions of NGC 4535 have no neutral gas (Cayatte et 
al. 1990). 

NGC 4536: This large barred galaxy displays two main promi- 
nent spiral arms and additional spiral-like structures. The only kine- 
matical perturbations are seen in the Z-shape of the velocities in the 
central parts and streaming motions along the arms. As in the CO 
data (Sofue et al. 2003a), a steep velocity gradient is observed in 
the Ha data in the galaxy core. 

NGC 4548: The observation has a very low initial S/N which 
causes the bin sizes to be large during the Voronoi tessellation. The 
disc rotation is nonetheless seen, with a velocity gradient of ~ 200 
km s^^, which result is in agreement with that deduced from long- 
slit spectroscopy (Rubin et al. 1999). 

NGC 4567: This galaxy is part of a pair with NGC 4568. The 
morphology and the velocity field do not show signs of disturbance. 

NGC 4568: Here again, no clear signs of velocity disturbance are 
observed, except perhaps in the external regions of the disc. Koop- 
mann & Kenney (2004) claim that this galaxy pair likely experi- 
ences RPS and a tidal interaction. It is however not obvious that a 
tidal interaction is (or has) occuring (occured) while regarding the 
very mildly perturbed kinematics of the pair. 

IC 3583: The galaxy presents irregular stellar and gaseous mor- 
phologies. Tschoke et al. (2001) revealed the presence of a faint, 
narrow spur of star formation regions to the South, pointing toward 
NGC 4569. This spur is also seen in the FP field. A weak rota- 
tion is evidenced in the perturbed velocity map (~ 100 km s^^ in 
projection). 

NGC 4569: The velocity field of the galaxy is perturbed, show- 
ing gas streamings along the arm-like structures in the Ha 
(severely) truncated disc. An off-plane string of Hll regions is ob- 
served to the West of the disc. It exhibits a velocity gradient that 
follows the rotation of the disc. This gaseous structure has already 
been seen in deep Ha image (Tschoke et al. 2001) and in Hi data 
(VoUmer et al. 2004). The giant diffuse gaseous region to the West 
of the disc (Tschoke et al. 2001) is not detected in the FP observa- 
tion, which is probably due to a too low S/N or to a likely cut-off 
in the transmission of the interference filter. Vollmer et al. (2004) 
proposed a ram pressure stripping model to explain the Hi asym- 
metry. 

NGC 4571: This almost face-on galaxy has a very patchy distri- 
bution of Hll regions. The interference filter unfortunately cuts part 
of the emission in the NE region. The 21 -cm data show an Hi dis- 
tribution well confined in the optical disc boundaries (van der Hulst 
et al. 1987). The Ha kinematics is in agreement with the Hi data 
and presents a slightly perturbed velocity field. 

NGC 4579: The FP data presents an Ha kinematics of NGC 
4579 (M58) which is perturbed by the streaming motions along the 
large-scale spiral arms and by a rather large velocity dispersion in 
the receding half (by comparison with the approaching side). The 
circumnuclear region is particularity interesting because it exhibits 
a bright "loop" structure within a radius of ~ 1.5 kpc (Keel 1983; 
Pogge 1989; Gonzalez Delgado & Perez 1996) and the nucleus is 
known to have an extremely broad Ha profile with a full-width at 
zero intensity of up to 18 000 km s~^ (Barth et al. 2001). The so- 
called "loop" is actually an asymmetric nuclear spiral which dis- 
plays a very large velocity gradient (~ 500 km s^^). The kine- 
matical major- axis of the nuclear spiral is almost perpendicular to 
that of the large-scale spiral structure. The high-resolution data of 
the molecular gas component showed CO emission concentrated 
within the core (Heifer et al. 2003; Sofue et al. 2003a) and the 



neutral hydrogen emission was lying in an annulus-like structure 
(Cayatte et al. 1990) which is actually the large scale spiral arms. 

NGC 4580: This Ha truncated disc displays little kinematical ir- 
regularities and a small velocity gradient. Though the Ha covers 
the innermost parts of the stellar distribution like in NGC 4405, a 
constant velocity seems to be reached along the major-axis, which 
was not seen in NGC 4405. The Ha morphology is in the form of a 
pseudo-ring. No ionized gas is detected in the spiral arms that can 
be seen in the broad-band images. 

NGC 4639: This barred galaxy exhibits a remarkable inner ring 
(Evans et al. 1996; Gonzalez Delgado et al. 1997). The outer 
gaseous distribution is furthermore reminiscent of an outer pseudo- 
ring. This outer structure is visible in the _B-band image but hardly 
in the NIR image, which principally shows the bar and the bulge. 
The velocity field is rather regular, consistent with a flat circular 
velocity pattern of the disc. It also shows streaming motions in the 
ring. A steep velocity rise is seen in the galaxy core, as was also 
observed in Rubin et al. (1999). 

NGC 4654: The stellar and gaseous morphologies are very asym- 
metric. The Ha velocity field is strongly perturbed. Alternating pat- 
terns of (decreasing and increasing) solid-body and constant veloc- 
ity are observed through the field, which is well illustrated by the 
PV diagram. In the approaching half, the transition region between 
the central solid body part and the outer decrease occuring at an 
offset of —1' leaves a remarkable imprint in the velocity field. In 
addition, streaming motions are evidenced along the spiral arms. 
A solid-body rotation is partly associated with the very thin and 
elongated Ha bar-like structure in the central parts. According to 
Vollmer (2003), a scenario mixing effects of a tidal interaction with 
NGC 4639 ~ 500 Myrs ago and RPS could be responsible of the 
perturbed stellar and gaseous distributions of NGC 4654. There is 
a large difference of distance between the two galaxies, ~ 22 Mpc 
and ~ 25 Mpc for NGC 4639 as estimated from the TuUy-Fisher 
relation (Yasuda et al. 1997) and the Cepheids method (Saha et al. 
1997) respectively, and ~ 14 Mpc for NGC 4654 (TuUy-Fisher re- 
lation, Yasuda et al. 1997). This inconsistency of distance points 
out that the Tully-Fisher distance of NGC 4654 must be highly er- 
roneous if the Vollmer's scenario is valid. 

NGC 4689: As for NGC 4571, the Ha emission is very patchy. 
The velocity field appears slightly perturbed, showing evident signs 
of streamings associated with the pseudo-spiral arms. The Ha kine- 
matics is in good agreement with the Hi (Cayatte et al. 1990) and 
CO (Sofue et al. 2003a) data. 

A2 Observational data 

Due to the restricted size allowed for 
the astro-ph files, only the example 
of NGC 4254 is presented here. All 
other 29 images can be previewed on 
http://www.astro.umontreal.ca/fantomni/virgoJ 
Full ps (36 Mb) and pdf 
(3Mb) resolution versions with 
all images are available on 
ftp://ftp.astro.umontreal.ca/outgoing/cheminl 
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Figure Al. Presentation of the Ha maps of the 30 Virgo galaxies. NGC 4254 (VCC 0307). Top : GOLDMine B-band (Left) and H-band (Right) images. 
Bottom : FP Ho integrated map (Left) and Ha velocity field (Right). All three images of integrated intensities are displayed using a square root stretch with 
increasing intensity levels from blue to red colours. In the Ha flux and velocity maps, a cross gives the position of the photometric centre. The velocity field 
is displayed using a linear scale and are in kms-i.ThePVdia gram is displayed using a logarithmic stretch. A solid line represents the derived rotation curve 
of the galaxy (see section|l}. 1' corresponds to 4.7 kpc. 



